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Abstract. Signals of bimodality have been investigated in experimental data of quasi-projectile decay 
produced in Au+Au collisions at 35 AMeV. This same data set was already shown to present several 
signals characteristic of a first order, liquid-gas-like phase transition. For the present analysis, events are 
sorted in bins of transverse energy of light charged particles emitted by the quasi-target source. A sudden 
change in the fragmentation pattern is observed from the distributions of the asymmetry of the two largest 
fragments, and the charge of the largest fragment. This latter distribution shows a bimodal behavior. The 
interpretation of this signal is discussed. 

PACS. 05.70.Fh Phase transitions: general studies - 25.70.-z Low and intermediate energy heavy-ion 
reactions 



1 Introduction 

The existence of different phases for infinite nuclear mat- 
ter is predicted by theoretical calculations since the early 
80 's [I]. Then, the possibility of observing a nuclear liquid- 
gas phase transition in the laboratory has been deduced 
from several experimental observations associated to the 
multi- fragmentation of finite nuclei. These observations in- 
dicate the occurrence of a state change in finite nuclei, 
which is interpreted to be the finite system counterpart of 
a phase transition [2]. 

Many of these signals are qualitative, and therefore 
cannot give information on the detailed trajectory (in terms 
of pressure, volume, temperature, isospin) followed by the 
system from one phase to the other. Other signals give in 
principle quantitative information, but can be distorted. 
Indeed the products of the reactions are detected asymp- 
totically and not at the production time, and therefore 
they need to be corrected for secondary decay [3]. These 
corrections are, at least partially, model dependent and 
induce systematic errors which are difficult to estimate 
quantitatively. 

To overcome these difficulties, it is important to per- 
form a systematic study of different phase transition sig- 
nals. The best would be to exploit new generation An ap- 
paratuses, in order to be able to investigate several signals 
at the same time, with the same experimental data sam- 
ples, and with a complete or quasi-complete detection [4]. 
Waiting for these new apparatuses, some of the signals in- 
dicating a phase transition have been obtained with mea- 
surements performed by the Multics [5] , together with the 
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Miniball [6 , multi-detectors. In the last few years we have 
investigated in detail the properties of quasi-projectiles 
detected in Au + Au reactions at 35 A.McV, with a fixed 
source charge, and at different excitation energies. The 
following signals have been obtained: 

1. the average size of the heaviest fragment (tentatively 
associated to the liquid part) decreases for increas- 
ing excitation energy of the nuclear system [7] with 
a power law distribution of exponent /3 ~ 0.31; 

2. temperature measurements result compatible [8,9 with 
a "plateau" in the caloric curve [TP] : 

3. critical exponents have been extracted [7j, close to the 
values expected within the liquid-gas universality class; 

4. the size distribution presents a scaling a la Fisher |llj : 

5. interaction energy fluctuations, corrected for side-fee- 
ding, were shown to overcome the statistical expec- 
tation in the canonical ensemble, corresponding to a 
negative branch of the microcanonical heat capacity 
for a system in thermodynamical equilibrium [12 . 

Some of these signals are coherent with the findings of 
other experimental collaborations with different data sets 
[13 lfT4lfT5l ll6j. In particular, the last two signals have been 
confirmed in central reaction measurements performed with 
the Multics [T7UT5] and with the Indra [9, 19] apparatuses. 
Some of these behaviors were also observed in other finite 
physical systems undergoing a transformation interpreted 
as a first order phase transition, namely in the melting of 
atomic clusters j201l2T| and in the fragmentation of hydro- 
gen clusters [2"2"] . 

Recently [23 , a new topological observable has been 
proposed to recognize first order phase transitions. When 
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a finite system undergoes such a transition, the most prob- 
able value of the order parameter changes discontinuously, 
while the associated distribution is bimodal close to the 
transition point, i.e. it shows two separate peaks, corre- 
sponding to the two different phases 24,25]. In the case 
of transitions with a finite latent heat, this behavior is 
in agreement with the Yang-Lee theorem for the distribu- 
tion of zeroes of the canonical partition sum in the com- 
plex temperature plane [26], and equivalent to the pres- 
ence of a curvature anomaly in the microcanonical entropy 

Since many different correlated observables can serve 
as order parameters in a finite system, the task is to choose 
an order parameter which can be accessible from the ex- 
perimental side [28]. This is the case for observables re- 
lated to the measured charges. The INDRA collabora- 
tion [29] has proposed as order parameter the variable 
Z S ym = zl+zt ' wnere Zi and Z2 are the charge of the 
largest and the second largest fragments detected, in each 
event, in the decay of an excited source. An indication 
of a bimodal distribution was obtained for this quantity. 
Signals of bimodality in different observables have been 
obtained in experiments with different projectile-target 
combinations, and in different energy ranges [13 ] f28 l l30] . 

In ref. [28] it has been pointed out that the variable 
lisym can present a spurious bimodality in small three- 
dimensional percolation lattices close to the percolation 
threshold. This behavior is due to finite size, and makes 
a bimodality in 7j sym an ambiguous signature of the tran- 
sition. On the other side, the size A\ or charge Z\ of the 
largest fragment have distributions which for any lattice 
size are consistent with the critical percolation behav- 
ior [3T]. These observables were then suggested as more 
apt to discriminate between a first order phase transition, 
a critical phenomenon, and a smooth cross-over. 

In this paper we investigate if these different signals 
are present in our data. We also discuss the relation of the 
bimodality signal with other phase transitions indicators 
obtained for our data. 



2 The experiment 

The measurements and the analysis have been extensively 
described elsewhere [7]. Here we recall that the measure- 
ments were performed at the K1200-NSCL Cyclotron of 
the Michigan State University. The Multics [5J and Mini- 
ball [6 j arrays were coupled to measure light charged par- 
ticles and fragments with a geometrical acceptance of the 
order of 87% of Att. The events have been recorded if at 
least two different modules have been fired. Similarly to 
Ref. [25] . the selection of the quasi- projectile (QP) source 
have been made by a shape analysis, keeping all the frag- 
ments with Z > 3. The fragments have been considered as 
belonging to the QP, if forward emitted in the ellipsoid ref- 
erence frame. The complete source has been obtained by 
doubling the forward emitted light particles in the back- 
ward direction, in order to minimize the contamination of 
light particles emitted by a possible mid-velocity source. 



At the end of this procedure, only events with total charge 
within 10% of the Au charge have been considered for the 
analysis, in order to study the decay of a well detected 
constant size source, in a wide range of excitation ener- 
gies. 
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Fig. 1. Charge density distribution of QP Au+Au events, as a 
function of the fragment velocity, along an axis parallel to the 
QP velocity. The continuous line (full points) represents the 
charge density for fragments accepted for the QP, the dashed 
line (open squares) for fragments of the QT, the dotted line 
(open circles) is a filtered simulation of a QT source symmetric 
to the experimentally detected QP. 



In order to visualize the source characteristics in the 
selected events, the fragment (Z> 3) charge density distri- 
bution [32] is shown in Fig. [1] as a function of the fragment 
velocity in the QP reference frame. The ensemble averaged 
charge density {pz(v par )) is defined as 

(pz{Vpar)) = ^ ^rj ) 

where pz(v par ) is the event-by-event distribution in 
the velocity v par for the collected charge fraction. This ob- 
servable represents the distribution of the collected charge 
bound in fragments along the direction of the QP veloc- 
ity. In Fig. [TJ the continuous line (full points) represents 
the charge density for fragments accepted for the QP, the 
dashed line (open squares) for fragments of the QT. This 
latter is consistent with the filtered simulation of a QT 
source symmetric to the QP (dotted line - open circles). 
QP and QT can be easily recognized, showing that the im- 
posed conditions are effective to select events where the 
contamination of a mid-velocity source is negligible (for 
more details see Refs. 0(9] )• 

The characteristics of these events have been examined 
by analyzing the isotropy of the fragment angular distri- 
bution in the quasi- projectile reference frame [T2"Ill7| . and 
by comparing the data to predictions of a statistical mul- 
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Fig. 2. Left panel: charge of the largest fragment Zi normal- 
ized to the size of the source Zo as a function of the QP calori- 
metric excitation energy. The symbols refer to event selected 
with (open circles) and without (full points) the constraint 
on the velocity of the largest fragment. The dashed line is a 
power law with exponent /3 = 0.31. The right panel shows the 
normalized partial energy fluctuations for QP events selected 
with (grey contours) and without (full points) the constraint 
on the velocity of the largest fragment. The estimation for the 
canonical heat capacity Ck is also shown |12] , 



tifragmentation model [38] . The general conclusion is that 
an important degree of equilibration appears to be reached 
by the excited quasi-projectile sources in the whole range 
of excitation energies. For more details, see Refs. [7] 9,12, 

EJ. 

In Refs. [TIE] a further condition to characterize the 
QP was added, i.e. that the velocity of the heaviest frag- 
ment is larger than 75% of the projectile velocity. This 
condition indeed has been replaced with the "complete- 
ness" of the event [281129] which reflects on the limitations 
on the parallel momentum, since it results less correlated 
to the variables we want to study. This does not affect 
the distribution shown in Fig. [TJ and all the signals of 
phase transition do not change significantly. To quantify 
this statement, we present in Fig. [5] the power-law be- 
havior of the average charge of the largest fragment, nor- 
malized to the charge of the source, as a function of the 
excitation energy, and the normalized partial-energy fluc- 
tuations, leading to the estimate of a negative branch for 
the microcanonical heat capacity [12]. The power law in 
the Zi distribution and the partial energy fluctuations are 
very scarcely affected by the condition on the velocity of 
largest fragment, apart from some very small variations 
at the higher energies [51112]. 

These observations mean that the detection appara- 
tus appears especially effective in the complete detection 
of purely binary collisions. The same is not true for all 
reaction mechanisms, which may need different detection 
systems to be addressed. Only a fraction of well detected 
peripheral collisions can be interpreted as the independent 
statistical decay of two isotropic sources [33] . For instance 
within the INDRA apparatus it has been pointed out that 
for 80 A.MeV Au+Au collisions, these events represent 
about the 20% of the total number of complete events [34] , 



and depend on the selection criteria adopted [35] . In our 
case the statistical events represent about the 30% of the 
measured events as reported in Ref. [36], and as can be 
inferred from Fig. 2b) of Ref. [7]. The different in the per- 
centage of statistical events could also be due to the differ- 
ent trigger conditions of two and four modules fired, used 
in our and Indra measurements, respectively. 



3 Signals of bimodality 

In the liquid-gas phase transition, the largest fragment 
detected in each event is a natural order parameter, be- 
cause of its correlation with the particle density in the 
grancanonical ensemble [231137] , The variable proposed in 
Ref. [291139] in turn is trivially correlated to the largest 
cluster size, and in addition brings further information 
on the global fragmentation pattern. This means that it 
should be possible to observe a bimodal distribution for 
the charge of the largest fragment or the asymmetry, if one 
considers a system close to the transition temperature. 

The global distributions of Z± and Z SJ/m are shown in 
Fig. [3] for all QP events selected as explained in section 
2. Because of the impact parameter geometry, this distri- 
bution is clearly dominated by peripheral collisions at low 
deposited energy, leading to a heavy Zj w 75 residue with 
a large asymmetry Z sym « 0.9. 




25 50 75 0.5 1 

Fig. 3. Distribution of the largest fragment charge Zi, (left 
panel) and of the asymmetry between the two largest fragments 
'Zsym (right panel) for the whole set of QP events. 



We can however also notice the presence of a large 
bump, corresponding to much lighter residues and much 
more symmetric fragmentation patterns. For this bump to 
be interpreted as an indication of bimodality, we should 
show that: 

— the two different decay patterns can be obtained in the 
de-excitation of the same source, 

— they correspond to the same temperature. 

Let us first concentrate on the source definition. As we 
have already stressed in the last section, we are consider- 
ing only events with a detected charge in the forward QP 
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hemisphere close to the original Au charge. This guaran- 
tees a good detection, but does not constrain the reaction 
mechanism or the number of sources, since the system is 
symmetric. Fig. Q] shows that the selected events are con- 
sistent with a purely binary kinematics, meaning that the 
bump at low charge shown by Figj3] cannot be ascribed 
to a reduced size of the excited source. However Fig. [1] 
is obtained with the whole set of events, which is largely 
dominated by peripheral collisions. We may then wonder 
whether a (small) contamination of central collisions, lead- 
ing to an important stopping in the center of mass, may 
be responsible of a decrease of the QP source size in the 
dissipative reactions corresponding to the low Z± bump. 

In Fig. Q] we plot the velocity in the laboratory frame of 
the QP source as a function of the excitation energy, with 
cuts of Zi > 50 (full points) and Zi < 50 (open points). 
The source velocity expected for a two body (QP-QT) 
kinematics, obtained via energy and momentum conser- 
vation in the hypothesis of an equal sharing of the excita- 
tion energy by the two collision partners, is shown by the 
dashed line [7] . We can see that for all calorimetrically re- 
constructed excitation energies, and for both Zi cuts, the 
observed behavior is compatible with purely binary col- 
lisions. This comparison shows that also the lighter Zi's, 
corresponding to constant size QP remnants, do not come 
from central reactions. 



sent two different phases, this means that the associated 
phase transition should have a non zero latent heat, as 
it is the case for regular liquid-gas. Therefore, the sort- 
ing variable should not impose a too strong constraint on 
the deposited energy, such that the two phases can be 
accessed in the same bin. In particular, in the liquid-gas 
phase transition, Zi is known to be bimodal in the canoni- 
cal ensemble which allows huge energy fluctuations, while 
no bimodality is observed in the microcanonical ensemble 
with fixed energy |23j . 

To search for a possible bimodal behavior, we should 
then in principle sort the data in temperature bins, i.e. in 
canonical ensembles. This is not possible experimentally, 
but we can choose a sorting variable as close as possible to 
a canonical temperature. Moreover, as suggested by pre- 
vious papers [291159] , the sorting observable should better 
not be auto-correlated with fragments and light particles 
emitted by the QP source. To fulfill these requirements, 
as in [2J5], the transverse energy Eti2 = Y2z Ezsin 2 (6z) 
of the light particles (Z < 2) emitted by the quasi-target 
source has been chosen, which is only loosely correlated to 
the QP observables. The QT had been much larger than 
the QP, this sorting could be considered as a canonical 
one. This is not our case, but it has still the advantage 
to allow for relatively large energy fluctuations, as it is 
needed to explore two phases that could be separated by 
a non zero latent heat. 




z (A.MeV) 



Fig. 4. Velocity of the QP source as a function of the excitation 
energy. Full and open points refer to Zi > 50 and < 50, respec- 
tively. The dashed line is obtained by energy and momentum 
conservation in a two body QP and QT kinematics [?]• 



This discussion implies that we can safely consider the 
data as characteristic of the de-excitation of constant size 
source in a wide range of excitation energies. 

Let us now come to the central question of data sort- 
ing. The global distributions of Fig.[5Jreflect the excitation 
energy deposit imposed by the dynamics of the entrance 
channel, and cannot be considered as belonging to a single 
statistical ensemble. If a sorting cannot be avoided, it is 
also clear that the shape of the distributions will depend 
on the sorting choice. 

The two de-excitation modes visible in Fig. [3] are asso- 
ciated to very different excitation energies. If they repre- 
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Fig. 5. Transverse energy distribution of light charged parti- 
cles emitted in the backward hemisphere in the c. m. reference 
frame. 



It has to be noted that efficiency limitations for frag- 
ments in the backward direction do not allow to have a 
total reconstruction of the quasi-target source (see Fig. [T|) . 
Therefore the selection on QT characteristics can only be 
performed on light particles; in particular the transverse 
energy Eti2 of these particles is detected with good effi- 
ciency The distribution of Eti2 for the selected events is 
shown in Fig. [5] 

This sorting can be assimilated to an impact parameter 
sorting. The width of the transverse energy bins, equally 
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spaced, was chosen of 40 MeV. Only the first six bins have INDRA collaboration on peripheral Xe + Sn and Au + 

sufficient statistics to be considered for the subsequent Au collisions [29] , 

analysis. 
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Fig. 6. Calorimetric e* distribution in the six bins of trans- 
verse energy Eti2. Eti2 increases going from left to right and 
from top to bottom panels. 
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Fig. 7. Charge of the heaviest fragment vs. the asymmetry 
Zsym of the two heaviest in different bins of transverse energy. 



The excitation energy constraint implied by this sort- 
ing is explored in Fig. [6l which shows the distribution of 
the calorimetric e* in the six transverse energy bins. We 
can notice from this figure that the variables Eti2 and e* 
are loosely correlated, and a relatively wide distribution 
of e* is obtained in most of the bins of transverse energy. 
It is well known [9] that the calorimetric measurement is 
not perfect, and the incomplete detection creates a spuri- 
ous width in the energy distribution. Since this spurious 
width never exceed 1 MeV per nucleon, it is clear from 
Fig[6] that the sorting in Eti2 bins cannot be considered 
as a microcanonical selection, where no bimodality would 
be expected. 

The charge of the heaviest fragment Zi is represented 
as a function of the asymmetry Z sym of the two heaviest in 
Fig. [7] We can see that the maximal probability does not 
monotonically change with the centrality selection. The 
most probable fragmentation pattern, characterized by a 
residue exhausting most of the available charge and an 
important asymmetry between the two largest fragments, 
abruptly changes between the fourth and the fifth bin with 
the apparition of a second peak. This peak represents mul- 
tifragmentation events, with the largest fragment compa- 
rable in size to the other emitted clusters. In the sixth bin 
this second peak tends to become more prominent, even 
if the situation would be more clear with higher statistics. 
Such a discontinuous behavior agrees with the expecta- 
tions from a phase transition, and with the findings of the 
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Fig. 8. Distribution of the Zi variable in the different trans- 
verse energy bins. 



Projecting the plots of Fig. onto the two axes, we 
can note that the best indication of bimodality appears 
on the Zi variable (see Fig. HJ, whereas the plot of the 
asymmetry Z sym does not show a clear bimodal behavior 
(see Fig. [9|. The largest fragment size distribution, peaked 
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around Z\ w 70 up to the fourth bin, shows a maximum 
around Z\ « 30 in the sixth bin, passing through a con- 
figuration (fifth bin) where a minimum in the probabil- 
ity appears to be associated to the intermediate patterns, 
even if the statistics should definitely be improved. This 
strongly suggests a first order phase transition [23] . 





Fig. 9. Distribution of the asymmetry Z ay m 
transverse energy bins. 



in the different 
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Fig. 10. Evolution with Eti2 of the most probable value of 
Zi (left panel) and Z sym (right panel). The lines are drawn 
to guide the eye. In this figure also the seventh bin of the 
transverse energy is shown, despite of the low statistics. 



A better understanding on the nature and order of 
the observed phase change can be achieved form Fig. [TTl 
which shows the distribution of the largest cluster charge 
and excitation energy in the Eti2 region (fourth + fifth 
bins) where the sudden change in the fragmentation pat- 
tern is observed. As we have already mentioned, a first or- 
der phase transition should imply a non zero latent heat, 
meaning that the two "phases" observed at the same "tem- 
perature" should be associated to different excitation en- 
ergies. 



The indication of bimodality is less clear in the 7* syrn 
variable, which shows a wide distribution in the last bins 
without any pronounced minimum ( see Fig. [5]). However 
the sudden change in the most probable fragmentation 
pattern, shown in the bidimensional correlations of Fig. [7] 
is clearly seen in both variables. The most probable value 
of Zi and 7i sym is shown in Fig [10] as a function of the 
transverse energy. Both variables show a sudden decrease 
from the evaporation dominated pattern, to the multifrag- 
mentation dominated one. This behavior is again in agree- 
ment with the findings of ref. [29] and consistent with the 
expectation from a phase transition. In principle a first 
order phase transition should be associated to a discon- 
tinuous jump [23] . while a continuous power law behavior 
would characterize a second order phase transition. This 
however would be true only if the sorting variable could 
be assimilated to a thermodynamical temperature. In the 
microcanonical ensemble, even a first order transition is 
associated to a continuous behavior of the order parame- 
ter. As a consequence, the power law behavior of the av- 
erage size of the largest cluster as a function of excitation 
energy (see Fig. [5J, can be observed both in the case of a 
critical behavior and in the coexistence zone of a first or- 
der phase transition [371140] . In the experimental case, the 
arbitrariness of the sorting and the absence of a physical 
external bath does not allow to draw definite conclusions. 




Z,/Z £ (A.MeV) 



Fig. 11. Distribution of Zi/Z (left) and E* (right) in the 
fourth + fifth bins of transverse energy. Thick lines: global dis- 
tributions. Thin lines: E* < 3.5A.MeV (left); Z x > 50 (right). 
Dashed lines: E* > 3.5 A. MeV (left); Zi < 50 (right). 



We can see in Fig.[TT]that indeed the two decay modes 
observed in the Z\ distribution correspond to different val- 
ues of the calorimctric excitation energy: the cut Zi=50 
that roughly identifies the two modes (see Fig. [3]) succeeds 
in splitting the energy distributions into two separate com- 
ponents (right part of Fig. [TT]) . even if the distance of the 
centroids is too small to produce a clear bimodality in the 
excitation energy distribution. The lower (higher) Z\ com- 
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ponent roughly corresponds to a deposited energy higher 
(lower) than 3.5 A.MeV (left part of Fig. [U]). This behav- 
ior is consistent with the expectations from a first order 
phase transition smoothed by finite size effects. 

If we interpret the two Z\ distributions as two coex- 
isting phases, it would be tempting to estimate the latent 
heat of the transition from the energy distance between 
the two peaks. The "liquid" peak points to an excitation 
energy E\ m 2 A.MeV which nicely agrees with the global 
Z\-E* correlation shown in the left part of Fig. [21 and 
with the indication of the fluctuation measurement shown 
in the right part of the same figure. The "vapor" contri- 
bution peaks at E% ~ 5 A.MeV, a lower value respect 
to the location of the second divergence in the fluctua- 
tion analysis. This discrepancy may be due to the intrinsic 
limitations of the transverse energy sorting, that does not 
allow sufficient energy fluctuations respect to a physical 
heat bath; it may also point to an incomplete exploration 
of the high energy phase space in our data sample, that 
cuts the distributions on the high energy side. 

To summarize, the results of Fig[Tl] tend to suggest 
that the observed sudden change from evaporation to mul- 
tifragmentation can be associated to first order phase tran- 
sition. Higher statistics samples obtained with collisions at 
higher beam energy could allow to be conclusive about the 
compatibility between fluctuations and bimodality |35] . In 
addition a detailed study of the convexity properties of the 
distributions is needed [HPS] , 



4 Conclusions 

In this paper we have presented a new analysis of the 
35 A.MeV quasi-projectile Au+Au data collected with 
the Multics-Miniball apparatus. The distributions of the 
largest cluster charge and of the charge asymmetry be- 
tween the two largest clusters detected in each event have 
been studied. These data allow to analyze, with a lim- 
ited statistics, the de-excitation of a constant size quasi- 
projectile source within a large range of dissipated energy. 
A clear transition from an evaporative to a multifragmen- 
tation pattern has been observed. The shape of the distri- 
butions have been studied to search for a possible bimodal 
behavior, which would allow to interpret this transition in 
the de-excitation mode as the finite system counterpart 
of a first order phase transition. The asymmetry distribu- 
tion does not present a clear structure, while the largest 
fragment charge appears bimodal. 

The same data have shown several different signals 
that coherently point to a first order liquid-gas-like phase 
transition. We recall here the determination of thermody- 
namically consistent critical exponents, both in a moment 
analysis [7] and in an analysis "a la Fisher" [TT], and the 
fluctuation peak in the partial energy distribution, with 
an absolute value of fluctuations consistent with the ex- 
istence of a negative branch for the microcanonical heat 
capacity [9TT2"]. 

Concerning bimodality, the search of a convexity in 
the distribution of the largest fragment emitted in each 



collisional event appears a very direct, and therefore in- 
teresting signature of a first order phase transition. The 
intrinsic weakness associated to this signal, namely the 
arbitrariness in the choice of the sorting variable and of 
the sorting bin width, will be the object of future investi- 
gations |41j . To confirm the significance of these findings, 
it will be important to verify that the dynamics of the 
entrance channel and the different bias associated to dif- 
ferent detectors do not influence the final results. To this 
aim, the bimodality distributions for systems produced 
with different entrance channels and detected with differ- 
ent apparatuses should be compared in details. Moreover, 
the stability of the signals as a function of the different 
ways of sorting the events should be checked. 

In particular, the qualitative trend of the data appears 
consistent with the recent findings of the INDRA collabo- 
ration [29] at higher incident energies, but the behavior of 
the asymmetry variable Z sym is not equivalent in the two 
data sets and deserves a further analysis. 
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